BEDT-TTF organic superconductors: the entangled role of phonons 
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We calculate the lattice phonons and the electron-phonon coupling of the organic superconductor 
K-(BEDT-TTF)2i3, reproducing all available experimental data connected to phonon dynamics. 
Low frequency intra-molecular vibrations are strongly mixed to lattice phonons. Both acoustic and 
optical phonons are appreciably coupled to electrons through the modulation of the hopping integrals 
(e-LP coupling). By comparing the results relevant to superconducting k- and /3*-(BEDT-TTF)2l3, 
we show that electron-phonon coupling is fundamental to the pairing mechanism. Both e-LP and 
electron- molecular vibration (e-MV) couplings are essential to reproduce the critical temperatures. 
The e-LP coupling is stronger, but e-MV is instrumental to increase the average phonon frequency. 



o 
o 

Of 

3 ; 



-a 

G 

o 
o 



> 



(N 
O 
<N 
O 

-i— > 
c3 



i 

C 

O 

o 



X 



Organic superconductors (oSC) have been studied for 
more than twenty years, but the nature of the pairing 
mechanism is still unknown, and remains one of the most 
intriguing problems in this class of materials Jf|. Like 
high T c superconductors, oSC are low-dimensional sys- 
tems characterized by strong electronic correlations. This 
fact and the proximity of the superconducting phase to 
magnetically or charge ordered states has brought about 
the suggestion of exotic pairing, mediated for instance 
by spin or charge fluctuations j|. However, direct ex- 
perimental evidence for such kind of purely electronic 
mechanisms has not been provided. On the other hand, 
the oSC phonon structure is very complex, preventing 
a clear assessment of the role of phonons in the pair- 
ing. Indeed, scattered experimental evidences of the in- 
volvement of intra-molecular phonons H , of optical inter- 
molecular phonons Q , and of acoustic phonons || have 
been presented. However, all these studies, focusing on a 
particular type of phonon, failed to offer a comprehensive 
and convincing picture of the phonon mediated mecha- 
nism in oSC. Aim of this letter is to provide a unified 
framework for the description of the phonon dynamics 
and of the electron-phonon coupling in oSC based on bis- 
ethylendithio-tetrathiafulvalene (BEDT-TTF). We shall 
show that all the above types of phonons are coupled to 
electrons, and cooperate in establishing electron pairing. 

In molecular crystals it is convenient to start by dis- 
tinguishing between intra- and inter-molecular vibra- 
tional degrees of freedom, and to consider their inter- 
actions at a later stage. Thus, in the framework of a 
tight-binding description of the electrons, intra-molecular 
vibrations modulate the on-site energy and yield the 
electron-molecular vibration ( e-MV) coupling, while lat- 
tice phonons, mainly inter-molecular in character, mod- 
ulate the hopping integrals (e-LP coupling). 

Molecular vibrations and e-MV coupling constants are 
transferable among different crystalline phases, and ei- 
ther computational or experimental methods can be used 
for their determination. Indeed, BEDT-TTF molecular 
normal modes and e-MV coupling constants are fairly 
well known B . The problem of lattice phonons is much 



harder to deal with, because the experimental character- 
ization of lattice phonons and of e-LP coupling constants 
is difficult and subject to considerable experimental un- 
certainty j|. Moreover, the phonon properties are not 
transferable among different crystals. To overcome these 
difficulties, we have adopted the "Quasi Harmonic Lat- 
tice Dynamics" (QHLD) computational method 0, ||, 
which is known to accurately predict the crystal structure 
and the phonon dynamics of complex molecular crystals 
as a function of pressure p, and temperature T. 

As extensively discussed elsewhere |), in QHLD 
the vibrational contribution to the crystal Gibbs en- 
ergy G(p, T) is approximated by the free energy of the 
phonons calculated in the harmonic approximation. The 
structure at a given p, T is determined self-consistently 
by minimizing G(p, T) with respect to lattice parameters, 
molecular positions and orientations. The method relies 
on an accurate inter-molecular potential model, which 
we have described as the sum of empirical atom-atom 
potentials, known to be reasonably transferable among 
crystals containing similar molecules. We have restricted 
our attention to the (BEDT-TTF) 2 l3 salts, which consti- 
tute an ideal benchmark for our approach. Our strategy 
has been to determine and to fix once and for all the 
needed C, S, H and I atom-atom parameters on the basis 
of the crystal structures of neutral BEDT-TTF and non- 
superconducting a-(BEDT-TTF) 2 I 3 salt g. We have 
then applied the method to the phonon structure and 
e-LP coupling strength of superconducting /3*-(BEDT- 
TTF)2l3 salt Jt]]. Here we illustrate the corresponding 
results for k-(BEDT-TTF) 2 I 3 . The comparison with the 
experimental data for the two superconducting /?*- and 
K-phases allows us to validate our approach and to assess 
the role of phonons in the superconductivity mechanism. 

Since both IJ and BEDT-TTF present low frequency 
intra-molecular modes falling in the region of the lattice 
phonons, the coupling between inter- and intra-molecular 
phonons may be important for the (BEDT-TTF)2l3 salts. 
To account for this effect, we use the exciton-like model 
described elsewhere 0. We stress that lattice phonons 
are actually mixed only with low-frequency vibrations 
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FIG. 1: Calculated (cv, full line) and experimental (c p , cir- 
cles, from Ref. [jl]) k-(BEDT-TTF) 2 I 3 specific heat. 

(below 300 cm -1 ), and these do not modulate on-site 
energies Therefore e-LP and e-MV couplings can be 
dealt with separately. 

At room temperature, k-(BEDT-TTF) 2 I3 crystallizes 
in the monoclinic system P2\/c, with four molecules per 
unit cell. Upon cooling below 150 K the crystal struc- 
ture changes to the P2 X system [§. The BEDT-TTF 
layer is in the 6c plane. The QHLD calculations repro- 
duce both crystal structures, which correspond to two 
distinct G(p, T) minima. The calculated lattice parame- 
ters are within 3% of the experimental ones. The results 
presented below all refer to the low temperature K-phase. 

A comparison with K-(BEDT-TTF)2l3 scanty infrared 
and Raman data does not offer a stringent test for the 
calculated phonon frequencies, whose number largely ex- 
ceeds that of the currently observed bands. We then 
compare the experimental specific heat c p to the cal- 
culated cy, which depends on the overall distribution of 
phonon frequencies uiqj-. cy (T) = fee ^ FChuo^/k-QT), 
where F(x) = x 2 e x (e x — l) -2 and the sum is extended to 
all phonon branches j and wavevectors q. As shown in 
Fig. [I], the agreement is outstanding. We stress that there 
are no adjustable parameters in our cy(T) calculation. 

To calculate the e-LP coupling strength we assume 
that the low frequency phonons, mainly intermolccular 
in character, couple to electrons only through the mod- 
ulation of hopping integrals: 

t KL =t° KL + J29(KL;qj) (1) 

3 

where t^L is the hopping integral between neighboring 
pairs KL of BEDT-TTF molecules, and Qqj are the di- 
mensionless normal coordinates. Finally, g(KL; qj) = 
•v/ h/ 2 ujqj {dtKL / dQ^j) defines the linear e-LP coupling 
constants. The t° KL and their dependence from the 
QHLD eigenvectors Qqj have been calculated by the ex- 
tended Hiickel (EH) method 0. The obtained t° values 
are in good agreement with those calculated by extended 
basis set density functional methods jllj]. The t° are 
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FIG. 2: Dispersion curves of k>(BEDT-TTF) 2 I 3 acoustic lat- 
tice phonons along the principal directions shown in the inset. 
Inset: Fermi surface calculated with the dimer model. 

used to describe the band structure in a tight binding 
approach, adopting the dimer model p2| . The result- 
ing conduction bands and Fermi surface (inset of Fig. |^) 
match closely those obtained considering all the four site 
orbitals within the unit cell . 

To express the e-LP coupling constants in reciprocal 
space, we have to introduce the dependence on both elec- 
tronic (k) and phononic (q) wavevectors. We have fol- 
lowed the approximation scheme adopted for /3*-(BEDT- 
TTF)2l3 0. The k dependence has been accounted for 
in terms of the dimer model: the coupling constants 
associated with the modulation of the intra-dimer hop- 
ping integral are assumed independent of k, whereas the 
inter-dimer g's are assumed to have the following depen- 
dence: g(k,k';j) = J2r9kl{<v) (e lk ' R - e lkR ), where 
q = k' k, and R represents the nearest neighbor lattice 
vectors. For what concerns the q dependence, the optical 
modes are taken independent of q, whereas for acoustic 
phonons we have simply assumed gKh{<^j) oc \J\ q | up 
to the average value calculated at the zone edges. 

Fig. |2| shows the frequency dispersion of the acous- 
tic modes along representative Brillouin zone directions, 
together with a pictorial representation of the coupling 
strength at the zone edges. The dots diameters are pro- 
portional to the squared coupling constant, which are 
much larger along the Y, Z directions {i.e. the &*, c* 
axes) than along the diagonal. Our results indicate that 
for k- (BEDT-TTF) 2I3 the coupling of acoustic phonons 
is very anisotropic. 

From the coupling constants in reciprocal space 
and the QHLD phonon density of states we calculate 
the Eliashberg coupling function a 2 {u>)F{uj) — 

^(^)E J (lg(k,k';j)| 2 <5(^ - Wqi ))FS, where N{E F ) is 
the unit cell density of states/spin at the Fermi level, and 
( ) fs indicate the average over the Fermi surface. Fig. || 
compares k- and /5*-(BEDT-TTF) 2 I 3 Eliashberg func- 
tions. Both curves have been smoothed by a Gaussian 
of constant width to ease the comparison. Although sig- 
nificant differences exist, also due to the higher number 
of phonons for the K-phase unit cell, the most significant 
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FIG. 3: Calculated Eliashberg function of k-(BEDT-TTF) 2 I 3 
(lower curve) and /T-(BEDT-TTF) 2 I 3 (upper curve, offset for 
clarity). The shaded areas at the bottom indicate the posi- 
tion of the most strongly coupled phonons, as detected from 
Raman on others K-phase BEDT-TTF salts (from Ref. Eh. 




FIG. 4: Eigenvector of k-(BEDT-TTF) 2 I 3 zone edge (Y) 
acoustic phonon at 21 cm -1 . 



peaks are correlated, as indicated by the dotted lines. 
In both salts the lowest frequency peak results from the 
overlap of acoustic and lowest frequency optical modes. 

Comparison of the computed Eliashberg function with 
available experimental data provides a test for the e-LP 
constants and QHLD eigenvectors. For the (3 phase, 
the calculation compares quite favorably j?j with the 
Eliashberg function directly derived from point con- 
tact measurements. This kind of measurement is not 
available for k-(BEDT-TTF) 2 I 3 , but we can make ref- 
erence to the phonon self-energy effects measured by 
Raman on k-(BEDT-TTF) 2 Cu(NCS) 2 and k-(BEDT- 
TTF) 2 Cu[N(CN) 2 ]Br [|. A quantitative comparison is 
not possible, due to the different counter-ions. Never- 
theless, Fig. H shows that the frequency regions (shaded 
areas at the bottom) of the few strongly coupled optical 
phonons |GJ, indeed correspond to the peaks of k-(BEDT- 
TTF) 2 l3 Eliahsberg function. We finally mention that in- 
elastic neutron scattering on k-(BEDT-TTF) 2 Cu(NCS) 2 
H evidences also the involvement of acoustic phonons, 
in agreement with our findings. 



Although further measurements, like Raman spectra 
with excitation in the infrared, are highly desirable to 
provide further tests, we believe that the above compar- 
ison already gives confidence to the QHLD eigenvectors, 
and as an example we show in Fig. Ej the eigenvector 
corresponding to K- (BEDT-TTF) 2 I 3 zone edge acoustic 
phonon at 21 cm -1 . As it occurs in several other cases, 
the coupling is induced by the relative displacement of 
BEDT-TTF along the long molecular axis, i.e. roughly 
normal to the BEDT-TTF layers. Moreover, a normal 
mode coupled to electron is in general a mixture of both 
BEDT-TTF and I3 vibrations. This mode mixing indi- 
cates that also the counter-ion vibrations may play an 
indirect role in the coupling. 

We now proceed to estimate the critical temperatures. 
From Fig. || Eliashberg function we derive the dimension- 
less coupling constant A LP , and the logarithmic average 
frequency |||: 
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The A value depends on N(Ep). With our tight binding 
EH calculation we estimate N(Ep) — 5.4 and 3.8 spin 
states/(eV unit cell), so that A LP = 0.59 and 0.65 for k- 
and /3*-(BEDT-TTF) 2 I 3 , respectively. The correspond- 
ing ajj^ values are 27 and 40 cm -1 . When these values 
are inserted into the Allen-McMillan equation 113]: 



l c = — exp 
1.2 v 



1.04(1 + A) 



A- M *(1 + 0.62A) 



(4) 



assuming the standard 0.1 value for the Coulomb pseu- 
dopotential one obtains T c — 0.8 and 1.7 K for k- 
and /3*-(BEDT-TTF) 2 I 3 , respectively, to be compared 
with the experimental values of 3.4 and 8.1 K |Q, [LtJ . 
Of course, the T c values depends critically on N(Ep), a 
parameter difficult to evaluate [Q. However, even by 
taking the maximum current estimates, derived from ex- 
periment and therefore enhanced by many body electron- 
electron and electron-phonon effects [[l4j], the T c values 
(2.7 and 3.6 K) remain significantly below the experi- 
ment. Therefore e-LP coupling alone cannot reasonably 
account for the T c of the two salts. 

However, one must not forget high frequency e-MV 
coupled phonons j^, ^| . By including these phonons into 
the calculation of A and uji n , the above scenario changes 
substantially. We underline that e-MV coupling gives an 
additive contribution to the total A, A = A LP + A MV , 
whereas the contribution to wi n is multiplicative: Wi n = 



;, LP c^ v , where 



,MV 



exp 



2N(E F ) 
NX 



with 



gi the e-MV coupling constants and N the number of 
molecules per unit cell 15 . It turns out that by in- 



cluding the e-MV coupling, A increases by only ~ 25 %, 
but uj\ n reaches 60 and 94 cm -1 for k- and /?*- (BEDT- 
TTF) 2 I 3 , respectively. By using these new values into 
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Eq. (3), keeping \i* = 0.1, we find that our tight bind- 
ing EH N(Ep) estimates, which are within the range of 
the currently accepted values H], account for the ob- 
served critical temperatures. In particular, precise T c 
matching requires N(Ep) = 5.2 and 3.9 spin states/(eV 
unit cell), corresponding to a total A of 0.74 and 0.91 
for k- and /3*-(BEDT-TTF) 2 I 3 , respectively. We finally 
remark that only by using the new uj\ n , which does not 
depend on N(Ef), the empirical relation [|l6] connect- 
ing this parameter to the specific heat jump AC p /"/T c 
yields for k-(BEDT-TTF) 2 I 3 a jump in nice agreement 
with experiment (1.7 vs 1.6) 

In our approach electron-phonon coupling provides a 
consistent and plausible explanation of many experimen- 
tal findings related to the superconducting properties of 
K- and /?*-phase BEDT-TTF salts. Besides the just men- 
tioned T c values and specific heat jump, the scattered ex- 
perimental evidences of the involvement of acoustic and 
optical lattice phonons have been rationalized. The cou- 
pling of acoustic phonons in the K-phase turns out to be 
strongly anisotropic, likely yielding an anisotropic gap 
that might confuse the node search |l7|) . In addition, we 
have found that the most strongly e-LP coupled BEDT- 
TTF modes imply motions approximately normal to the 
conducting planes, with counter-ions vibrations mixed to 
BEDT-TTF ones (cf. Fig. |). This finding might well 
account for the "interlayer effects" detected in thermal 
expansion measurements jL8|. Experimental evidences 
in favor of the involvement of e-MV coupled phonons are 
mainly related to isotopic effects || |l!|. In particular, 
the universal inverse deuterium isotope effects (0.25 K) is 
likely due to a small increase in u>i a . The e-MV coupling 
constants of deuterated BEDT-TTF are not known, but 
we have verified that very small and plausible changes of 
the coupling constants with respect to those of pristine 
BEDT-TTF || may indeed lead to the required increase 
in T c , despite the large downwards shift of several fre- 
quencies upon deuteration. 

The principal conclusion of our paper is that phonons 
are the main responsible for the coupling mechanism in 
BEDT-TTF based oSC. The role of phonons is com- 
plex, or "entangled" , in the sense that different kinds 
of phonons are all important, but contribute in different 
ways. The acoustic and optical lattice phonons, modu- 
lating t, give the main weight to A, but intra-molecular 
phonons modulating on-site energies are essential in in- 
creasing the average wi n value. Although here we vindi- 
cate the role of phonons, we also stress that traditional 
phonon mediated pairing alone cannot give reason of 
all the experimental findings related to superconducting 
properties of BEDT-TTF salts. The proximity of antifcr- 
romagnetic or charge ordered states in many oSC clearly 
indicate that electron-electron interactions are important 
, and these cannot be properly accounted for in our 
effective single particle approach. In particular, just con- 



sider the above mentioned k-(BEDT-TTF) 2 Cu(NCS) 2 
and K-(BEDT-TTF) 2 Cu[N(CN) 2 ]Br salts. Their iu ln is 
similar to k-(BEDT-TTF) 2 I 3 , yet their T c 's are 9.5 and 
11.5 K, respectively ||(|. We think it is unlikely that the 
corresponding increase in A with respect to k- (BEDT- 
TTF)^ can be attributed to a significant change in the 
e-LP coupling strength due to the change in counter- 
ions. Rather, we believe that electron-electron interac- 
tions and/or the proximity of a phase boundary may 
lead to an enhancement of the effective electron-phonon 
coupling. In other words, we need to analyze the role 
of electron-phonon coupling in low-dimensional strongly 
correlated materials, still a poorly understood topic fel] . 
Our work represents the necessary prerequisite for such 
an analysis in BEDT-TTF-based oSC. 
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